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- Abstract 

"We present the analysis of an XMM-Newton observation of the M 17 nebula. The X-ray point source population consists of massive 
O-type stars and a population of probable low-mass pre-main sequence stars. CEN la,b and 01352, the X-ray brightest 0-type 
stars in M 17, display hard spectra (kT of 3.8 and 2.6 keV) consistent with a colliding wind origin in binary/multiple systems. We 
show that the strong interstellar reddening towards the O-type stars of M 17 yields huge uncertainties on their Lx/L^oi values. The 
low-mass pre-main sequence stars exhibit hard spectra resulting from a combination of high plasma temperatures and very large 
interstellar absorption. We find evidence for considerable long term (months to years) variability of these sources. M 17 is one 
"of the few star formation complexes in our Galaxy producing diffiise X-ray emission. We analyze the spectrum of this emission 
■and compare it with previous studies. Finally, we discuss the Optical Monitor UV data obtained simultaneously with the X-ray 
'images. We find very little correspondence between the UV and X-ray sources, indicating that the majority of the UV sources are 
.foreground stars, whilst the bulk of the X-ray sources are deeply embedded in the M 17 complex. 
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'1. Introduction 
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' M 17 is one of the most luminous H n regions of our Galaxy. 
.This nebula hosts the very young open cluster NGC6618 
which probably reflects the second stage i n a still on-go i ng se 



quential (massive) star formation process (IHanson et all 11997 



Povich et al. 



2009'). make M 17 a privileged 

.target for the observational study of massive star formation. In- 
deed, M 17 harbors several cocoon stars, which are still embed- 



ON 



ded in the relics of their cra dles of c i rcumstellar material (e .g 



Chini&Kriigel. 1985; Hanson et all 119971; IChini et all BOOOc 



Kassis et alll2002h . Several interesting objects are found along 



"the interface between the H ii region and the molecular cloud 
.M 17SW, supporting the idea that they were formed through 
"the interaction between the H n region and the molecular cloud. 
. ^ .Recent studies reported evidence for still accreting massive or 
intermediate-mass protostars for a couple of objects which were 
H .identified as silhouette disks sometimes associated with Ht jets 



(e.g. Niirnberg eret all l2007t iNielbock et all 120071 120081) . 



The distance of the M 17 complex remains somewhat un- 
certain. From multi- color photometry and spectroscopy, 
iHoffmeis ter et al.1 (12008^ inferred a distance of 2. 1 + 0.2 kpc and 
a ratio of total to selective extinction of Ry = 3.9. These au- 
thors argued that previous photometric and spectrophotometric 
distance estimates, around 1.6 kpc, were probably biased by ig- 



" Based on observations collected with XMM-Newton, an ESA Science Mis- 
sion with instruments and contributions directly funded by ESA Member States 
and the USA (NASA). 
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noring the multipli city of O-type st ars in NGC6618 (see also 



the discussion in Povich et al. [ l2009l) 



The total extinction towards M 17 consists of a foreground 
extinction, a local dust component inside M 17 and, in some 
cases, absorption by circumstellar disks or dust. This results 
in a very la rge optical extinctio n with a very complex spatial 



distribution dPovich et al. [12009') 



The stellar population of M 17 features a number of very 
young, hot and m assive objec t s. Th e main exciting source 
of M 17 is CEN 1 JChiniet all Il980l) . also known as Klein- 
mann 's Anonymous Sta^ CEN 1 is a Trapezium-like system 
consisting of two spectroscopic binaries (CEN la and lb) that 
form a visual pair with an angular separation of 1.8arcsec. 
Each of these spectrosco pic binaries is of spectral type 04 V 



)pic 

B 



(iHoflrmeister et all 12008114 Actually, multipl city might be 
quite common among the O-type stars of NGC6618. Indeed, 
CEN 18 (06-8V) was also found to display double absorp- 
tion lines, and CEN 37 (03-6 V), being over-luminous for its 
spectral t ype, was suggested to be another binary candidate 
(iHoffmeister et all 12008.) . VLA radio observations at 3.6 and 
6 cm wavelength revealed several compact radio sources, two 
of them being asso ciated with the co mponents of the visual bi- 
nary CEN 1 (Rodriguez et al. [ l2009h . CEN la is detected only 
at the shortest wavelength. Both radio sources, especially the 
one associated with CEN la, have flux densities that exceed the 
expected level of free-free emission from their stellar winds. 



' T hroughout this paper, we adopt the naming c onvention of IChini et alj 
i\98(l CEN designations) and lOgura & Ishid j <1976l 01 designations). 

^The physical association of the two bi naries is still contro versial since they 
are subject to different optical extinctions iHoffmeister et ali[2008l) . 
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Rodriguez et al. I (l2Q09h accordingly suggest that the observed 
radio emission could, at least partially, stem from the syn- 
chrotron process in a colliding wind region of the binary sys- 
tem. 

Although stellar winds account only for a very small fraction 
(a few percent at most) of the overall energy input of early-type 
stars into H ii regions, their mechanical power determines the 
dynamics of a nebula and produces so-called wind-blown bub- 
bles or superbubbles (in the case of clusters of early-type stars). 
The superbubble blown by the stellar winds of the OB stars in 
M 17 is highly asymmetric due to its interaction with the neigh- 
boring molecular cloud M 17SW. Soft X-ray emission spatially 
coincident with the superbubble of M 17 was detected with pre - 
vious X-r ay observations using ROSAT dPunne etall l2003b 



Chandra jTownslev et al.L 12003a b) and Suzaku dHvodo et al 



I2OO8.) . Beside this diffuse emission, these X-ray observations 
also revealed a number of point sources. For instance, a 40 ksec 
Chandra-AClS observation led to the detection of 886 point- 



like sources in the field of M 17 (Broos et all 120071) . 



In the present paper we discuss the results of an XMM- 
Newton observation of M 17 with the goal to quantify the prop- 
erties of the brightest X-ray sources in the field. Section|2]pro- 
vides an overview of our observations. In Sect.[3] we analyze 
the X-ray sources found in the field of M 17. We cross-correlate 
their position with various catalogs to identify counterparts at 
other wavelengths. We then analyze the spectra of twenty X- 
ray bright point like sources as well as of the diffuse X-ray 
emission. Finally, in Sect.|4] we summarize the most important 
results of our study and present our conclusions. 




Figure 1: Energy-coded thi'ee-color X-ray image of Ml 7 built from the EPIC- 
MOS data. North is up and east is to the left. The field of view of the EPIC 
instruments is centered on CEN 1 and has a radius of 15 arcmin. The red, green, 
and blue colors correspond to the soft, medium, and hard energy bands used 
throughout this paper (see Sect.|2)- The individual images were exposure cor- 
rected, before they were combined. Note the reddish emission lane that extends 
to the east-south-east of the central cluster 




Figure 2: Energy-coded three color UV image of M 17 as built from the UVW\ 
(red), UVM2 (green) and VVW2 (blue) OM data. The OM field of view covers 
a region of 17 X 17 arcmin^, centered on CEN 1. Note the diffuse emission that 
is seen eastward of the center The brightest point-like sources are affected by 
some artifacts (see Sect.[2)- 



2. Observations 



M17 
satellite 



was observed 
( Jansen et al.[ 



for 



2001) 



(JD 2 452 7 10.170 - 24 5 2710 .520) 
de tectors dTurner et al. 



dStriider et al 



30 ksec with the XMM-Newton 
on 3 November 2003 
The two EPIC-MOS 
and the EPIC-pn camera 



2001) 



200 lb were all used in full frame mode with the 



medium filter to reject optical light. 

The raw data were processed with the SAS software (version 
10.0). Our observation was affected by a high radiation back- 
ground. As the background showed no isolated flares, but rather 
a general high level aff'ecting the entire duration of the observa- 
tion, we did not define good time intervals and kept the entire 
observation instead. 

Images were extracted for soft (0.5 - l.OkeV), medium (1.0 
-2.0keV), and hard (2.0- 10. OkeV) energy bands. An energy- 
coded three-color image of the combined field of view of the 
two EPIC-MOS instruments is shown in Fig. [T] In this image, 
we have emphasized the soft band to highlight the weak dif- 
fuse emission that extends to the east-south-east of the central 
cluster. 



The Optical Monitor (OM, iMason et al.L 1200 Ih instrument 



aboard XMM-Newton observed M 17 through the three UV fil- 
ters UVWl, UVM2 and UVW2 with exposure times of 4, 4 
and 8 ksec respectively. The OM was operated in 'full frame 
low-resolution' mode covering a 17 x 17 arcmin^ square field 
of view centered on the same coordinates as the X-ray cam- 
eras. The UVWl, UVM2 and UVW2 filters have bandwidths 
of 2450 - 3200, 2050 - 2450, and 1800 - 2250 A, respectively 



(iMason etaUl2001h . 
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Figure 3 : Color-magnitude (left) and color-color (right) diagrams of the near-lR counterparts of the X-ray sources detected with EPIC in the M 1 7 field of view. Only 
2MASS data with reliable photometry are plotted. The known O-stars and the bright X-ray sources are objects that are studied spectroscopicall v in Sect. | 3.2l The 
slope of the reddening relation in the color-color diagram, indicated by the dashed lines, is taken to be E{J - H)/E(H — K) = 2in agree ment with lHoffmeister et alj 
(2008) and Straizvs & Laugalvs (2008). The solid fine yields the locus of the un-reddened main-sequence according to lTokunagal i200dl) and lMart ins & Plez (2006) 
and shifted by a distance modul us of 1 1.61 in the co lor-magnitude diagram. The dotted straight line in the color-color diagram indicates the locus of the un-reddened 
classical T-Tauri stars following lMever et al.l ^9^). 



The OM data were processed with the relevant command 
chain of the SAS software. The data from the three filters were 
used to build a color-coded image (Fig.|2]i. The OM images ex- 
hibit some artifacts such as smoke-ring ghosts caused by light 
reflected off a chamfer in the detector window and modulo-8 



pattern around bright sources dMason et al.L 1200 Ih . Neverthe- 
less, Fig. |2] reveals also some genuine diffuse emission in the 
UVWl filter. A total of 614 point sources are detected in the 
OM images (a Table with all the detections will be made avail- 
able electronically), 128 of them being detected in all three 
passbands. 

The OM detector has a read-out time of 1 1 ms and sources 
with a count rate of 10 coun ts s~' typical l y suff er from coinci- 



dence losses of about 10% ( Mason et all 1200 lb . These losses 



are partially recovered during the data processing to provide 
'correcte d' count rat es, which are then converted into AB mag- 
nitudes (iOkel Il974l) . With the filters used during our obser- 
vation, corrected count rates of lOcountss ' correspond to 
UVWl, UVM2 and UVWl AB-magnitudes of 16.066, 14.912 
and 14.072, respectively. The reconstructed photometry of 
sources much brighter than this should be considered with cau- 
tion. 



3. X-ray sources in M 17 

Figure[T] shows that the X-ray emission of M 17 is dominated 
by a few bright point sources along with a number of fainter 
point sources and some diffuse emission. The SAS detection 
routines produced a list of 99 sources in the combined EPIC 



images with a significance threshol(|^ of 10. 

All these sources were inspected by eye and 16 of them were 
rejected as spurious. This inspection revealed another 9 sources 
that were clearly seen on the image, but were not found by the 
detection algorithm. The reasons the detection routines did a 
rather poor job on these data are essentially the diffuse back- 
ground as well as the generally high level of particle back- 
ground. We ended up with a final list of 92 sources (a table 
with these detections will be made available electronically). 
This is far less tha n the number of sources (886) detected by 



Broos et al.l ( 120071) with Chandra. The main reasons for this 



situation are the wider point spread function of XMM-Newton 
compared to Chandra and the high background level of our ob- 
servation. From a histogram of the number of sources versus 
count rate, we conclude that our detections are complete down 
to rates of 0.004 and 0.002 cts s ' for the EPIC-pn and the indi- 
vidual EPIC-MOS cameras, respectively. 

3.1. Cross-correlation with other catalogs 

To start, we have cross-correlated the lis t of X-ray s ources 



with the 2MASS po int sources catalog dCutri etall 12003 



Skrutskie et alj, 20061) . Using the technique of Jeffries et al 



[si 

(11997), the optimal cross-correlation radius, which offers the 
best compromise between obtaining the maximum number of 
true identifications and avoiding contamination by spurious co- 
incidences, was established to be 2.5 arcsec. For this correla- 
tion radius, we expect to find fewer than 7% of spurious cross- 
identifications among the 67 positive matches listed. 



^^This implies a probability of < e that a random Poissonian fluctuation 
could have caused the observed source counts inside the detection cell. 
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Table 1 : Results of the cross-coiTelation between the list of X-ray sources detected with our EPIC data of M 17 and the list of OB stars from lPovich et al.l <2009l) . 
Spectral types with a colon are uncertain, and could actually be earlier than quoted here. 



XID 


cc (J 2000.0) 


5 (J 2000.0) 


Name 


Spectral type 




t/yw2 


UVM2 


UVWl 


# 


hh:mm:ss 


deg : arcmin : arc sec 






mag 








38 


18:20:25.9 


-16:08:33 


CEN18 


06-8V 


7.6 






18.85 + 0.05 


45 


18:20:27.5 


-16:13:32 


OI345 


06V 




16.12 + 0.02 


16.05 + 0.01 


13.94 + 0.01 


55 


18:20:29.8 


-16:10:45 


CENla,b 


04V + 04V 








19.73 + 0.11 


68 


18:20:34.5 


-16:10:12 


CEN2 


05V 


5.2 


15.95 + 0.01 


15.83 + 0.01 


13.77 + 0.01 


71 


18:20:35.4 


-16:10:49 


CEN3 


09V: 


3.7 


13.59 + 0.01 


13.51+0.01 


13.05 + 0.01 


72 


18:20:35.9 


-16:15:43 


OI352 


08V: 


7.0 


20.17 + 0.23 


20.19 + 0.19 


17.65 + 0.02 


88 


18:21:02.2 


-16:01:01 


BD-16°4826 


05 


3.9 









The J versus J -H color-magnitude and J -H versus H - K 
color-color diagrams of the 2MASS counterparts of the sources 
detected with EPIC are shown in Fig.[3 The main-sequ ence 



Tokuna ga (2000) for 



relation in these figures is taken from 
stars later than spectral type O, and from iMartins & Plez 
(20Q3) for O-type stars. As can be seen on these diagrams, 
there is a wide range of reddening values, which hampers a 
direct comparison with theoretical isochrones. We note that 
about 1/6 of the X-ray sources have a near-lR counterpart that 
falls in the region occupied by classical T-Tauri stars (cTTs) 
which exhibi t IR-ex cess emission from a circumstellar disk 
(IMever et all Il997h. This result is in excellent agreement 



with the finding of iBroos et al.l (120071) who noted that whilst 
most of the Chandra sources occupied the near-lR color space 
associated with reddened weak-line TTauri stars, about one 
sixth of their sources fell in a region of the near-lR color-color 
diagram where classical TTauri and Herbig AeBe stars are 
expected. This contrasts with the situation in some other 
young open cl usters (e.g. CygOB2, 'Rauw'201l'; NGC6231, 
Sana et alj2007; and NGC 6383. Rauw et al...2003l) . where only 
a very small fraction of the counterparts of faint X-ray sources 
provide evidence for near-lR excesses. 

We have cross-correlated our list of sources with the list of 



OB stars provided by iPovich et al.l (l2009h . 16 known O-type 
stars fall within the field of view covered by the EPIC instru- 
ments. Within a correlation radius of 2.5 arcsec, we found 
7 positive correlations (see Table[T]i. Except for CEN 3 and 
Ol 352, all of them have a spectral type earlier than 07. The 
non-detections mostly refer to stars later than 07. There are a 
few exceptions though. CEN 37 (03-6 V) is located very close 
to the pair CEN la,b and could contribute to the X-ray emission 
associated with the latter, without being resolved as an individ- 
ual source. CEN 43 (03-5 V) falls close (6.9 arcsec) to one of 
the brightest X-ray sources (# 64) in our observation, but their 
association is formally not significant. Whilst we cannot rule 
out that CEN 43 slightly contributes to the emission of source 
#64, it must be stre ssed though th at this star is very heavily red- 
dened (Ay = 12.3. IPovichetall [2009) and could thus escape 
detection simply as a result of the interstellar absorption. 

In addition, this cross-correlation identifies the two most 



heavily reddened sources in Fig. [3] as being two candidate 
young stellar objects (YSOs) previously seen with Chandra: 
our X-ray sources 17 a nd 75 being respect ively associated with 
objects X4 and X14 in lPovich etd] ( |2009|) . 



Correlation of our X-ray sources with the catalog of OM 
point sources yielded a rather low number of matches: only 1 1 
of the 76 EPIC sources that fall into the OM field of view were 
found to have an OM counterpart. Six of these objects are OB 
stars (see Table[T]). The poor match between the X-ray sources 
and the UV detections is yet another manifestation of the huge 
absorption towards M 17. 

This brings up the question about the nature of the OM 
sources. Cross-correlation of the OM sources with the 2MASS 
point source catalog yields an optimal correlation radius in the 
range 2 - 2.5 arcsec. Adopting a radius of 2.5 arcsec, 5 10 of the 
614 OM point sources have a near-lR counterpart. The corre- 
sponding J -H versus H-K color-color diagram as well as the 
near-lR and UV color-magnitude diagram (see Fig.|4| clearly 
indicate that the majority of these objects are less-reddened, 
likely foreground late-type stars, unrelated to M 17. 

Comparing the list of OM sources with the catalogs of OB 
stars and candidate young stellar objects of Povich et al. ( 20091) 
yields 8 positive matches for OB stars and one YSO within a 
radius of 2.5 arcsec. Two of the OB stars seen with the OM are 
not detected in X-rays. These are CEN 25 and 31 of spectral 
type 07-8 V and 09.5 V respectively. 

Finally, we have cross-correlated the list of our X-ray sources 
with th e X-ray sources de tected by the Chandra-AClS obser- 



vation dBroos et al.L 1200 7). Adopting a correlation radius of 
2.5 arcsec, we find 61 positive correlations. Actually, some of 
our brightest EPIC sources do not have an ACIS counterpart 
(within the correlation radius) and that some of the prominent 
ACIS sources are missing in our list of detections. We have 
then used these results to compare the ACIS and EPIC count 
ratefl Figure|5]reveals a highly scattered relation, that deviates 



"^For the latter we adopted the MOS instruments as they have fewer gaps 
between CCDs than the pn detector, thus allowing a reliable determination of 
the count rates for a larger number of sources. 
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Figure 4: Left: UV color-magnitude diagram of tlie OM sources detected in the UVWi and UVM2 filters. Middle and right: color-color and color-magnitude 
diagrams of the near-IR counterparts of the UV sources detected with the OM in the M 17 field of view. Only 2MASS data with reliable photometry are plotted. 
The slope of the reddening relation and the main-sequence relations are built in the same way as in Fig.[3] 



significantly from a simple proportionality. One can consider 
several explanations for the large scatter: 

• there could be a bias in the EPIC count rates due to source 
confusion. Indeed, the density of X-ray sources as seen 
with Chandra is quite high, especially in the crowded 
core of the NGC6618 cluster. Due to the coarser PSF 
of XMM-Newton (as compared to Chandra) several Chan- 
dra sources could be confused and counted as a single 
source by XMM-Newton. To test this scenario, we have 
repeated the cross-correlation with a radius of lOarcsec 
and we have replotted only those EPIC sources with a sin- 
gle ACIS counterpart within lOarcsec (see Fig.|5]l. This 
treatment should eliminate the bulk of the sources that are 
due to confusion, but we see little effect on the distribution 
of the objects in Fig.|5] Therefore, confusion alone cannot 
account for the scatter in the count rate relation. 

• the brightest sources could suffer from pile-up. Because 
of Chandra's very narrow PSF, this problem would mainly 
impact the ACIS-I count rates. However, for all but the two 
brightest sources in Fig.|5] pile-up should be of a few per- 
cent at most. Indeed, for the brightest and second brightest 
objects in this figure, we estimate pile-up levels of about 
12%, and 5% respectively. For sources with less than 
0.01 counts"', the pile-up level should be less than 1.2%. 

• there could be systematics due to the instruments. Whilst 
one cannot rule out some residual problems in the cross- 
calibration of the EPIC and ACIS-I instruments, it has to 
be noted that a major effort has been put into this task over 
the last decade. All major issues have been solved and 
remaining problems are mostly small. Of course, some 
of the scatter in the count rates relation could reflect the 
different sensitivities of the EPIC-MOS and ACIS-I detec- 
tors. However, this should be a rather small effect. 

• there could be intrinsic time variability of the faint sources 
on long time scales. The Chandra and XMM-Newton 



observations were taken twenty months apart. It seems 
quite plausible that the X-ray brightness of weak sources, 
which could be associated with low-mass pre-main se- 
quence stars, changes with time. In view of Fig.|5] and 
the above discussion, this appears to be most likely ex- 
planation of the scatter and deviation from a simple pro- 
portionality. A spectacula r illustration of thi s scenario is 
provided by source 396 of 'B roos et all feOO?*). This is the 



second brightest X-ray source in the Chandra observation, 
but is not detected in our observatior||. Conversely, our 
fourth brightest source (XID # 64) has no equivalent in the 
ACIS-I data, the nearest possible ACIS counterparts being 
very faint and located at more than 4.5 arcsec. 



3.2. Spectral analyses 

The EPIC spectra of twenty point-like sources were extracted 
and f itted using the xspec software (version 12.6.0, Arnaud[ 
1996h . These include all the O-stars detected with EPIC as 
well as the brightest objects (with EPIC-MOS count rates of 
at least 0.01 counts"') in the field of view. The source regions 
were chosen in such a way to include the majority of the source 
photons, whilst simultaneously avoiding overlap with nearby 
sources in the most crowded regions of the field. The back- 
ground spectra were evaluated from nearby circular regions, 
free of any obvious point source. 

For all sources, we tested absorbed powe r-law models as wel l 
as absorbed thermal plasma apec models femith et al.Ll200lh . 
For the latter, we first assumed solar abundances. However, for 
the best exposed spectra, including the bright CEN 1 system, 
this resulted in a rather poor fit especially of the Fe K line. We 
then let the global metallicity vary during the fit and , in most 
cases, it converged to a value 0.2 - 0.4 times solar dMernier , 



""The closest EPIC source to this object is our source #42 at a separation of 
4.9 arcsec. 
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Figure 5: Comparison of the EPIC-MOS and AC IS-I count rates for the X-ray 
sources in common between ' Broos et"ai] ^2007*) and the present paper Open 
and filled symbols stand for EPIC sources which have a single ACIS counterpart 
within a coirelation radius of respectively 2.5 and 10 arcsec. 



2012h . However, given the limited quality of most X-ray spec- 
tra, a more robust approach would be to constrain the abun- 
dances independently and then fit the X-ray spectra using fixed 
abundances. 



Garcia-Roias et alJ (l2007h studied the chemical composition 



of the M 17 Hii region and determined abundances of nine key 
elements (including He, C, N, O, Ne, and Fe). Compared to 
the solar composition (lAsplund et al. [ l2009h . most elements are 
found to have abundances ~ 0.7 times solar There are a few 
notable exceptions to this rule. These are He, which has a solar 
abundance, C, which is found to be 2.4 times solar and, most 
of all, Fe, which is found to have an abundance of 0.026 solar. 
Garcia-Roias et al.l (,2007 ) show that their results are generally 



in good agreement with those from the literature, including for 
the iron abundance. The abundances of He and C have no di- 
rectly measurable impact on the EPIC spectra of our sources: 
helium has no spectral features inside the EPIC energy domain, 
whilst carbon has some lines that fall at the low energy edge 
of the EPIC sensitivity range, but are of little relevance here, 
because of the huge interstellar absorption towards M 17. How- 
ever, iron is of major importance for the EPIC spectra. Using 
such a low Fe abundance in our spectral models actually results 
in extremely poor fits. After some trials, we thus decided to fix 
all metal abundances to 0.7 times solar, except for He, C and Fe, 
which were fixed at 1.0, 2.4 and 0.4 times solar respectively. 

3.2.1. X-ray spectra of 0-type stars 

For the O-type stars listed in Table[T] we have constrained the 
value of the interstellar neutral hydrogen column density using 
the relation between Nh and E{B - V) derived bv lBohlin et al 



(Il978h . For this purpo se, we adopt the va lues of Ay quoted in 
Table[T]and taken from Povich et al.l ( 2009lR accounting fo r the 



peculiar reddening law Ry - 3.9 (iHoffmeister et al.L 120081) . For 



CEN la and CEN lb, there are two different values of Av of re- 



spectively 10.2 and 13.5 mag (IHoffmeister et al.L 120081) . In our 



fits of the X-ray spectra of O-stars, we allo w for some excess 



absorption due to the ionized stellar wind (iNaze et al.L l2004l) 



or to some circumstellar absorption. Therefore, we choose the 
lower of the two Ay values for constraining the Nh of the spec- 
trum of source # 55 (corresponding to CEN la -i- lb). The re 



is no e stimate of Ay for OI345 in the catalog of iPovich et al 
12009). However, the source is detected in all three OM filters 
and by comparison with CEN 2 (which has a similar spectral 
type), we estimate an extinction of Ay ^5.3 for OI 3450. 

All O-type stars detected with EPIC have X-ray spectra that 
are best fit with a thermal plasma model. The majority of the 
O-type stars in our sample have intrinsically rather soft spectra 
with kT between 0.3 and 0.6 keV. This is in line with the overall 
properties of O-type stars as derived from CCD X-ray spectra 
( Naze , 20091) and is consistent with the general picture of intrin- 
sic X-ray emission due to instabilities embedded in the stellar 
wind. There are two notable exceptions to this rule. These are 
CEN la,b and OI352. The spectra of these stars are shown in 
Fig.|6] They both display an Fe K line and require a second, 
hotter, temperature of 3.8 and 2.5 keV for CEN la,b and OI 352 
respectively. It is wort h comparing our be st-fit spectral param- 
eters to those found bv' Broos etal] (l2007b with the ACIS-I in- 
strument aboard Chandra. 

Concerning CEN la,b, we must recall that our spectrum in- 
cludes both CEN la and CEN lb and is likely also contaminated 
by CEN 37, although the latter was found to be 60 times fainter 
than the pair CEN la,b with Chandra. These three objects were 
resolved by Chand ra (sources 536, 543 and 574 respectively in 



Broos et al.L 120071) . Our EPIC spectra indicate a prominent Fe 
K line (much stronger than in the ACIS-I spectrcH) that con- 
strains the temperature of the hotter plasma component quite 
well. This temperature (3.8 keV) is significantly lower than the 
values of the hotter components found in the ACIS-I spectral 
fits, which were 15 keV for CEN la and 10.4keV for CEN lb 
(IBroos et al.L l2007h . CEN 37 was rep orted to have an in trin- 
sically soft spectrum (kT = 0.6 keV, iBroos et al.l l2007h and 
should thus not contribute to t he hard emis sion. Temperatures 
above lOkeV as advocated bv lBroos et al.l (i2007) are very dif- 



ficult to understand within the conventional scenarios (wind 
instabilities or magnetically channeled winds for the intrinsic 
emission, colliding winds in binary systems; see .Giidel & N aze. 
20091 for a review) for the X-ray emission of massive stars 



and would be reminiscent of so-called yCas-like variables. 
Our new result is much more in line with typical emission 
from a colliding wind binary system as would be expected 
from a Trapezium-like system such as CEN la,b, displaying 



We note that the se Ay values are larger than those quoted by 
IHoffmeister et al. |<2008l) for the stars in co mmon among our sample. 

'We note that lHoflfmeister et al] j2008l) quote Ay = 4.3 and 4.0 respectively 
for OI 345 and CEN 2. 

^We note that the ACIS-I spectrum of CEN lb is likely affected by pile-up. 
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Table 2: Results of the fits of the EPIC spectra of the 0-type stars of M17 with models of the kind wabs*wind*vapec and wabs*wind*vapec (2T) . The interstellar 
neutral hydrogen column densities were frozen at the values given in column 2. The absorption by the ionized wind is accounted for using the model of lNaze et alj 
(2004). The abundances of the hot plasma were taken to be 1.0, 2.4, 0.4 and 0.7 times solar for He, C, Fe and all other metals, respectively (see text). A two 
temperature plasma model is only required for CEN la,b and 01352. In these cases, column six quotes the ratio of the emission measures of the hotter and cooler 
components. Column 7 provides the reduced x' of the fit as well as the number of degrees of freedom (in brackets). The observed and absorption-corrected X-ray 
fluxes are evaluated in the 0.5 - lOkeV energy range. 



Star Nh logA^wind kT^ m EM2/EM1 xi (d-o.f.) fx /T^^ 

10^^ cm"^ cm"^ keV keV 10"'"' ergcm"^ s"' 10"'"' ergcm"^ s"' 



CEN 18 


1.13 


22.1!*'"^ 


0.3+"^ 








1.13 (27) 


0.11 


1.25 


01 345 


0.79 


< 21.6 


0.6!"| 








0.91 (29) 


0.51 


3.22 


CEN la,b 


1.52 


21 3+°"' 


61+"°^ 


3.8^ 


.2 
.2 


0.70 + 0.06 


1.13(455) 


42.30 


129.60 


CEN 2 


0.77 


21 6+°"2 


30+"°^ 






0.91 (74) 


0.81 


9.28 


CEN 3 


0.55 


21 5+°-^ 


30+"°'' 








1.43 (49) 


0.28 


2.09 


01 352 


1.04 


21 70+"-^ 


0.6+-4 


2.6^ 


.3 
.3 


1.15 + 0.27 


1.05 (234) 


11.06 


22.87 


BD-16°4826 


0.58 


91 7+0.4 










1.21 (29) 


0.64 


4.28 



CENla.b 01352 




0.5 1 5 10 0.5 1 5 10 

Energy (keV) Energy (keV) 

Figure 6: EPIC spectra of CEN la,b (left) and 01352 (right) along with our best fit models (see Table|2). The lower panels illustrate the residuals in the form of 
contributions to the overall 
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non-thermal radio emission (iRodrfguez et al. ■ l2009h . The ob- 
served X-ray flux of CEN la,b in die XMM-Newton observa- 
tion (4.23 X 10^'^ergcm"^ s"') exceeds the combined value of 
CEN la and CEN lb in the Chandra observation by about 37%. 
This difference could result from the combination of various ef- 
fects, including pile-up of the ACIS-I spectra, some differences 
in the adopted energy ranges (0.5 - 8.0 keV for Chandra versus 
0.5 - lOkeV used here), as well as genuine orbital variability, 
which is a common feature among colliding wind massive bi- 
naries. 

OI352 is the second brightest source in our data. Its 
spectrum is dominated by a hard emission component (kT = 
2.6 keV). IBroos et al.l (l2007h reported a temperature of 4.4 keV 
for the hotter plasma component. Again, our temperature 
is lower, but the difference is less dramatic than for CEN 1 . 
There is however a rather large discrepancy in the observed 
fluxes between Chandra {fx 



-10 



4.44 X 10^'^'ergcm^^s"') 



and XMM-Newton {fx = 11.06 x 10"' ' ergcm"^ s"')- Such 
a difference could hardly be attributed to instrumental effects 
and would thus likely indicate genuine flux variability. The 
high plasma temperature, the high X-ray luminosity (for an 08 
star) and the possible variability all point at this system being a 
colliding wind binary candidate. 

Using the interstellar absorption corrected X-ray fluxes and 
an estimate of the bolometric fluxes of the O-type stars, we 
could in principle infer the value of log (/™^''V/boi) of the O- 
type stars in M 1 7 and compar e it to the typical value of -6.45 + 
0.5 1 inferred by iNaze (l2009h . In the present case, the main 



uncertainties on log(/™^''V/boi) stem from the highly uncer- 
tain bolometric luminosities. This is illustrated in Fig.|7] where 
we have adopted two different approaches for estimating /boi. 
First we use the observed V magnitudes, the Ay values adopted 
for the evaluation of N h, and the bolometric corrections from 
Martins & Pled (120061) . This yields an extremely large value for 
the bolometric flux of the CEN 1 a,b system which would then 
be X-ray under- luminous. Whilst this method is independent of 
the actual distance of the stars, it is extremely sensitive to any 
uncertainties in Ay. Alternativel y, we have adopte d the typical 
bolometric luminosities from Martins et al.l (120051) . assuming a 
distance modulus of 1 1 .61 and accounting for t he multiplicity 
of CEN la, lb and 18 (iHoffmeister et all l2008b as well as for 
the suspected multiplicity of OI 352. This method is hence sen- 
sitive to the uncertainty on the distance modulus as well as to 
the assumptions on multiplicity. Figure|7] illustrates the huge 
impact of the uncertainty on the bolometric flux in the case of 
CEN 1 a,b which would be highly X-ray over-luminous within 
the second scenario. Nevertheless, we can conclude that on av- 
erage and within the uncertainties described above, the X-ray 
luminositi es of the O-s tars in M 17 follow the standard scaling 
relation of Naze (l2009h quite well. 



3.2.2. Other point sources 

In addition to the O-stars, we have also analyzed the spectra 
of 13 relatively bright point sources in the field of M 17. Nine 
of them have a near-IR counterpart (see Table|3]l. All spectra 
were fitted both with thermal plasma models, adopting the same 
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Figure 7: Absorption corrected X-ray flux of the O-type stars versus bolomet- 
lic flux. The filled symbols yield the bolometric fluxes derived from the V 
magnitudes and the Ay values, whilst the open symbols yield the bolomet- 
ric fluxes obtained assuming the stars are typical main-sequence O-stars at 
a distance of 2.1 kpc. See Sect. l3.2.l1 for details. T he solid line yields the 
log(/™''''V/boi) = -6.45 relation from lNazj )2009l) . The dashed lines cor- 
respond to the dispersion about this relation. 



abundances as for the O-stars, and with non-thermal power-law 
models. The latter would be relevant for background AGN, un- 
related to M 17. 

The results of these fits are displayed in Table[3] The ther- 
mal fits indicate very high temperatures, exceeding 2.4 keV, for 
all sources. Furthermore, with one exception, the best-fit col- 
umn densities exceed 0.5 x 10^^ cm"^ and reach 3.4 x 10^^ cm"^ 
for source #33. This corresponds to E{B - V) - 5.86 o r 
Av = 22.8, adopting the Ry value of ' Hoffmeister et al.l (l2008h . 



Hence, it seems unlikely that our sample is significantly pol- 
luted by foreground stars and most of these sources probably 
either belong to M 17 or are background objects. Assuming that 
these sources are located at the distance of M 17 and that the 
thermal plasma models are adequate, the absorption-corrected 
fluxes correspond to X-ray luminosities in the range 4.6 x lO"" 
to 5.5 X 10^^ erg s"' . The upper limit of this range corresponds 
to the highly absorbed source #33 and is thus more uncer- 
tain because of the correction of the interstellar absorption. If 
we discard this object, the highest luminosity still amounts to 
2.2 X lO^^ei-gg-i 

The high temperatures, the fact that most sources are located 
in M 17 or behind, the range of luminosities as well as the loca- 
tion of the near-IR counterparts in the color-magnitude diagram 
(Fig. [3]) suggest that these sources could be pre-main sequence 
(PMS) stars. PMS stars reach their highest X-ray luminosities 
during flares and we have thus inspected the X-ray light curves 
of all the sources of our spectroscopic sample. Whilst most of 
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Table 3: Results of the fits of the EPIC spectra of non-O-star point-like sources in the field of M 17 with models of the kind wabs*vapec (upper part of the table) 
and wabs*power (lower part). In the thermal plasma models, the abundances of the hot plasma were taken to be 1.0, 2.4, 0.4 and 0.7 times solar for He, C, Fe and 
all other metals, respectively (see text). Column 4 provides the reduced of the fit as well as the number of degrees of freedom (in brackets). The observed and 
absorption-corrected X-ray fluxes are evaluated in the 0.5 - 10 keV energy range. 
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the sources do not show a significant, coherent variability be- 
tween the different EPIC cameras, source #37 undergoes a flar- 
ing event towards the end of the observation, although our data 
only cover the onset of the flare, with the count rate increasing 
by a factor ~ 25. 

For about half of the sources in our sample, a slightly bet- 
ter fit quality is obtained with an absorbed power-law model. 
These fits are listed in the lower part of Table[3] When fitted 
with thermal plasma models, these sources indicate rather high 
best-fit temperatures: in some cases, only a lower limit to kT 
could be obtained. 

One source that is of special interest is #52. It is the 
third brightest non-O-star source with an X-ray luminosity of 
2.2 X lO-'^ergs"'. This object has two Chandra counterparts, 
CEN 98 (photometi-ic spectral type B4 V. Chini et al., 1980) an d 
CEN97 (photometric spectral type Bl V, IChini et all Il980l) . 
both located at 1.68 arcsec from the XMM source. There is 
no UV counterpart detected with the Optical Monitor, which 
is consistent with a deeply embedded source, as suggested 
by the lar ge Nh and the Ay of 10.6 and 11.0 ma g quoted by 



Chini et al. (1980) for CEN 97 and 98 respectively. iBroos et al 
(l2007h derived observed X-ray fluxes of 5.0 x 10"* and 6.3 x 



10"'^ ergcm"^ s"', which combined together are slightly less 
than our value (1.83 x lO"'-' erg cm"^ s"'). This source there- 
fore apparently corresponds to a pair of early-mid B-type stars 
with an exceptionally large X-ray luminosity. Indeed, typical 



X-ray luminosities of early B main seq uence stars are ab out a 
factor ten below the result found here dNaze et al.l 1201 lb and 
many of the mid-late B stars detected in X-rays are likely bi- 
nari es with a low-mass PMS companion being the X-ray emit- 



ter (lEvans et al.L 1201 ll) . It would be highly interesting to re- 



investigate the spectral types of these stars and to look for mul- 
tiplicity or other peculiarities. 

3.3. The diffuse emission 

As stated in Sect.[T] soft diffuse X-ray emission from the 
M 17 superbubble w as detec ted with p revious X-ray obs erva- 
tories ( Dunne et"all J2003 ; Townslev et a l. , 2003 a; Hvod o et al.i 
l2008h . This emission is strongest around the core of the ion- 
izing cluster, NGC6618, and further extends toward s the east- 
south-e ast du^ection as an elongated structure (Townslev et al.l 
2003b[ see also Fig.[B . 



iDunne et"aL (2003) noted a huge discrepancy between the 
observed level of X-ray emission and the two orders of mag- 
nitude higher level expected for a m odel of a wind-blow n su- 
perbubble including heat conduction dWeaver et al. , 1977 , see 
also Anorve-Zeferino et al.,2009). They suggested that this dis- 
crepancy stems from the effect of heat conduction and found 
a better agreement for models where heat conduction is sup- 
pressed (probably by magnetic fields), but dynamical mixing 
of cold interstellar material with the hot gas of the bubble is 
included. 
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Figure 8: MOS 1 and MOS 2 spectra of the diffuse X-ray emission in M 17. 



Because of their superb angular resolution, the Chandra 
data rule out the possibiUty that a substantial fraction (more 
than about 8%) of the diffuse emission coul d be due to un 



resolve d late-type pre-main sequence stars jTownslev et al 



2003bj). These data yielded an (absorption corrected) X-ray lu- 
minosity of - 3.4 X 10^'' ergs"' in t he 0.5 - 2keV band 
assuming a distance of 1.6 kpc ( Town slev et al. , l2003bl) . A 
very similar luminosity (Lx = 3.5 x 10^^ erg s~' in the 1.0 - 



S.OkeV energy band) was obtained bv lHvodo et al.l (l2008h from 
a Suzaku observation. 



Reves-Iturbide et al. (l2009l) presented 3-D hydrodynamical 



simulations of the M 17 superbubble and its X-ray emission. 
They successfully reproduced the observed champagne flow 
morphology and luminosity of the X-ray emission. These 
authors demonstrate the importance of the inhomogeneous 
interstellar medium and of the interaction of the cluster wind 
with the dense molecular cloud for correctly reproducing the 
observed X-ray luminosity. 

We have first extracted the spectrum of the diffuse emis- 
sion around the core of NGC 6618, excluding all detected point 
sources. This spectrum can be fitted with a two-te mperature 
thermal plasma model with kT of 0.8 and ~ 4.0 keV ( Mernieii 
2012!) . The higher temperature is most probably an indication 
of residual contamination by a population of unresolved point 
sources associated with pre-main sequence stars. 

As a next step, we have then extracted the diffuse X-ray 
emission within a rectangular box centered on (18:20:52.1, 
-16:12:09), with a size of 8 x 2.2arcmin and incUned by 20° 
with respect to the west-east direction, excluding two point 
sources. This b ox is slightly smaller th an the 8 x 3.5arcmin^ 
region used by iTownslev et al.l (l2003bh and avoids the cen- 
tral part of NGC 66 18 whilst extending somewhat farther 
away from the cluster The background was extracted over 
a circular area of radius 0.9arcmin centered on (18:20:33.1, 
-16:13:44). The resulting EPIC-MOS spectra are illustrated 



in Fig. [8] They are qualitatively similar to the spectrum pre- 
sented by iTownslev et al. ( 2003bh . The latter authors used a 
two-temperature model with one component (kT = 0.13keV) 
with solar composition and another one (kT - 0.6 keV) with 
non-solar abundances. The Suzaku observation was fitted with 
a single temperature ( kT - 0.25 keV) pl asma model with clearly 
sub-solar metallicity dHvodo et al. . 2008) . 

Here, we have adopted a different approach, fixing the 
plasma abundances of all components to the same mix as used 
for the point sources in M 17. The best-fit (xl - 0.89 for 
169 d.o.f.) parameters are Nh = 0.76+ |^ x lO^^cm ^ kTi = 
O.ll^-^^keV, kT2 = 0.22+12 kg V, eMz/EM, = 0.20 and fx = 
5.76 X 10"'^ ergcm"^ s"'. The observed flux is about half the 
value found with Chandra. This difference most likely reflects 
the somewhat different extraction regions used in the present 
study and in the analysis of the Chandra data. Our study thus 
reveals that the spectral parameters derived from the spectrum 
of the diffuse emission are quite sensitive to the treatment of 
several aspects, such as the plasma abundances, the definition 
of the extraction regions of the source and the background, and 
the possible contamination by unresolved point sources. 



4. Summary and conclusions 

In this paper, we have analyzed X-ray and UV data obtained 
during a 30ksec XMM-Newton observation of M 17. 

The brightest X-ray point sources are associated with seven 
O-type stars or systems of such stars. The X-ray images are 
dominated by the Trapezium-like system CEN la,b ((04 V -H 
O) -H (04V H- O)) and the 08V star 01352. In both cases 
we have shown that the EPIC spectra yield plasma tempera- 
tures that are si g nifican tly lower than found in a previous study 
by Broos et al.l (l2007h and are in very good agreement with 
a colliding winds scenario. An optical or near-IR spectro- 
scopic monitoring of these stars, although challenging because 
of the large interstellar absorption and unknown orbital peri- 
ods, would be extremely useful to establish their multiplicity 
and provide their orbital properties. 

There is a very large uncertainty on the Lx/^boi ratios of 
the O-type stars in M 17. This uncertainty mainly stems from 
the impact of the huge (and sometimes poorly known) redden- 
ing on the estimates of the bolometric luminosities. Whilst the 
Lx/Lboi ratios agree on average with the canonical relation es- 
tablished for a large sample of O-stars, the large uncertainties 
in the present case prevent us from a more detailed comparison 
that might reveal subtle differences and environmental effects. 

The majority of the secondary X-ray sources are likely low- 
mass pre-main sequence stars, with about 1/6 of the objects 
showing some evidence of near-IR excesses attributable to cir- 
cumstellar material. The near-IR counterparts of our X-ray 
sources span a wide range in reddening, which makes a quanti- 
tative comparison of their colors and magnitudes with theoret- 
ical evolutionary tracks very difficult. Out of the 13 brightest 
X-ray sources, not associated with an O-star, only one displays 
a flare in our 30ksec observation. This can be compared to the 
results of the 40 ksec Chandra observation: iBroos et alj (120071) 
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reported significant variability for 39 out of 886 sources, with a 
wide range of flare morphologies. Beside this short-term (flar- 
ing) activity, there must also be a substantial long-term (of or- 
der months or years) variability of the X-ray sources of M 17. 
Indeed, the relation between the ACIS-I count rates, recorded 
in March 2002, and the EPIC-MOS rates, measured twenty 
months later, exhibits a very important scatter which likely re- 
flects this variability. 

There is little correspondence between the 92 X-ray sources 
and the 614 UV detections. Most of the UV sources are in fact 
foreground objects, unrelated to M 17. The UV data also reveal 
some nebular emission eastwards of the central cluster 

Throughout this work, we have fitted the X-ray spectra of 
point and diffuse sources with plasma models accounting for 
the speci fic metal abun dances of the M 17 nebula as determined 
by Gai'cia-Roias et al ] tool ), except for the iron abundance. 
In deed, the latter is clearly h igher in our spectra than reported 
by Garcfa -Roias et al. I (l2007h . This is an interesting puzzle and, 
although an in-depth analysis of this effect is beyond the scope 
of the present paper, we suggest that this discrepancy might be 
related to ionization effects. The presence of a very hot plasma 
(as revealed by the diffuse X-ray emission) might impact on the 
ionization stages of iron and could bias the abundances inferred 
from optical diagnostics ([Feii] and [Feiii] lines). 
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